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Abstract. 

Based on the mathematically well defined Pade Theory, a theoretically safe new procedure for 
the extraction of the pole mass and width of resonances is proposed. In particular, thanks to the 
Montessus de Ballore's theorem we are able to unfold the Second Riemann sheet of an amplitude to 
search the position of the resonant pole in the complex plane. The method is systematic and provides 
QhI a model-independent treatment of the prediction and the corresponding errors of the approximation. 
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The non-perturbative regime of QCD is characterized by the presence of physical res- 
onances, complex poles of the amplitude in the transferred energy at higher (instead of 
the physical one) complex Riemann sheets. From the experimental point of view, one can 
obtain information about the spectral function of the amplitude through the Minkowsky 
region (q^ > 0) and also about its low energy region through the experimental data on 
I the Euclidean region (q^ < 0). In reference [1], the particular case of the nn vector form 

factor (VFF) was analyzed with the main purpose of studying its low energy behavior 
using the available Euclidean data. In particular, the first and the second derivatives of 
the VFF were determined at = with Pade Approximants (PA) centered at the origin 
trough a fit procedure to that data [1]. In such a way, the vector quadratic radius (r^)^ 
and the curvature Cy were extracted from the fit and, as a consequence, a value for the 
low-energy constant L9 = (6.84±0.07) ■ 10^-^ was obtained, [1, 2]. 

Despite the nice convergence and the systematical treatment of the errors, this proce- 
dure does not allow us to obtain properties of the amplitude above the threshold, such as 
in the case of the nn vector form factor, the p -meson pole position. The reason is simple: 
the convergence of a sequence of Pade Approximants centered at the origin of energies 
(q^ = 0) is limited by the presence of the n — K production brunch cut, [2] . The PA 
sequence converge everywhere except on the cut. Still, the mathematical Pade Theory 



allow us to produce a model independent determination of resonance poles when certain 
conditions are fulfilled. The most important one is to center our Pade approximant se- 
quence above the branch-cut singularity (beyond the first production threshold) instead 
of at origin of energies {q^ — 0). This small modification also provides the opportunity 
to use Minkowskian data in our study instead of the Euclidean one. The relevance of this 
model-independent method to extract resonance poles is clear since does not depend on 
a particular lagrangian realization or modelization on how to extrapolate from the data 
on the real energy axis into the complex plane. 

Although we apply this method in the particular case of a physical amplitude to extract 
the position of a resonance pole, is clear that it can be applied in a broader number 
of cases since only relies on a mathematical theory and not on a particular physical 
situation. We illustrate that method using an example where the properties above appear 
naturally. 

Imagine a function F{x) analytic in a disk B§{xq). Then, the Taylor expansion 
Fi?^) = 'Ln=o^n{x — xq)'^ Converges to F{x) in B§{xo) for AT — >■ oo, with derivatives 
a„ = F^"\xo)/n\. In that situation, one usually uses experimental data to extract the 
derivatives of F{x) using polynomial fits at higher and higher order N. Since the exper- 
imental data have errors, one normally finds that polynomials with order higher than a 
certain A^* do not produce new information, with the new coefficients of order A^* + 1 
been compatible with zero. In that situation, one stops the fit procedure to order A^* and 
believe that is the best one. The scenario changes, however, when the function F{x) 
is not analytic anymore, for example when has inside the disk B§{xo) a single pole at 
X — Xp. In this case, the Taylor series does not converge any more, and we need a different 
procedure to extract information about the function. 

The Montessus de Ballore's theorem states that the sequence of one-pole Pade Ap- 
proximants around xq, 

<(x,xo) = i' «*(x-^o)'+ Yz^fe^ , (1) 

converges to in any compact subset of the disk excluding the pole Xp, i.e, 

lim/f (jc,Jco) = F(jc). (2) 

As an extra consequence of this theorem, one finds that the Pade Approximant pole 
xpA = xo + converges to Xp iox N ^ °° as it must be as stated by Eq. (2). Since 
experiments provide us with values of Fj at different Xj instead of the derivatives of our 
function, we can use the rational functions as fitting functions to the data in a similar 
way as in Ref. [1]. In this way, as A^ grows P^ gives us an estimation of the series of 
derivatives and the Xp pole position. 

Usually, as we have already said, Pade Approximants are constructed around the low- 
energy point Xq = Q where x is the total energy squared. For a physical amplitude, 
the function F{x) (without a right-hand cut) is analytic from x = — oo up to the first 
production threshold Xfh and within the disk 5;c^^(0). In the %% vector form factor case, 
the threshold is found to be at Xj/, = 4m|, where niji is the mass of the pion. Experiments 
then may provide with data F^^p {x) at jc < and use them to extract the derivatives of the 



vector form factor at the origin [ 1 ] . One may also have Minkowskian F^^p [x + «0+ ) data 
at jc > which cannot be used to extract form factor properties with Fade Approximants 
centered at xq = due to the presence of the essential singularity at x = Xth- However, 
one can still use Fade Approximants in a safe way by using the Montessus' theorem 
and center the approximants at xq + /O"*" over the brunch cut between the first and the 
second production threshold, i.e., between Xth <xq < Xtu- In the nn vector form factor 
that would correspond to the range between pion production threshold and the kaon one, 
Xth = andff/, = 4mj. (assuming small multipion channels). In such a way we unfold 
the Second Riemann sheet due to the analytical extension of the function F(x) from the 
first Riemann sheet at x + iO^ into the second one. 

In the case of resonant amplitudes, a single pole appears in the neighborhood of the 
real x axis in the second Riemann sheet which can be related to a hadronic state, a res- 
onance. Once we have unfold the Second Riemann sheet by locating our approximants 
over the brunch cut, the application of the Montessus' theorem in the disk of conver- 
gence (which is the region defined between the thresholds) is straightforward and allows 
us to locate the position of the resonant pole if lies inside that region. If that is the case, 
our Fade approximant sequence determine systematically its position in a model 
independent way. 

In the next section we present the details of this procedure by a particular example 
and in the third section we apply all this technology to the real case of the ;r;r vector 
form factor ALEFH data to extract the mass and the width of the p -meson. 



TESTING THE METHOD WITH A MODEL 

To illustrate the possibilities of our method we consider a p-like model of the KK 
vector form factor, with a single pole in the second Riemann sheet at = (0.77 — 

/0.15/2)^GeV^ and a logarithmic branch cut (starting at = for simplicity). The 
model: 

has two parameters, M and F, tuned to produce the pole exactly at = q^. 

There are two different ways to explore the method with this model. The first one 
consists on using the derivatives of the model to construct a Pf Fade Approximants 
sequence and extract from them the convergence sequence for the position of the pole. 
We call this method the genuine method. The second method consists on simulating 
a physical situation generating a series of data points with zero error, which would 
represent an ideal experimental situation where all the uncertainty would be theoretical. 
We fit then that set of data for the modulus and the phase-shift of F{q^) with the 
corresponding modulus and phase-shift of each Fade Approximant and we extract the 
complex parameter a„ for each Pf . The Fade Approximant pole ^jj^ — {^fit ~ i^fit/'^)^ 
is found to converge very quickly to the real q^p = {Mp — iTp/2) ^ of the model as — > oo. 
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FIGURE 1. Rate of convergence of the sequence of genuine Pade Approximants corresponding to 
the first case studied for a p-like model (see the text for details). 



In the first case, using the genuine method, one first expand the function F{q^) at 
q-^ — q^, F{q^) = ao(^^ — ^o)'^ + ai(^^ — ^q)^ + ff{{q^ — q^Y), and with this expansion 
one construct the simplest Pade Approximant/'P((7^;^Q) = "''\_ . Then, one looks 

for the pole of this P^, q^ = "I" compares it with the pole position of the 

function F{q^). In a second step, one expands F{q^) at q^ = q-Q up to one order more, 
i.Q, F {q^) = aQ{q^ - q^)*^ + a\{q^ - q^y + a2{q^ - qQ^ + G'{{q^ - q^Y), and constructs 

the approximant pI (q^) = — ^ "21-2' — extract its pole q^^ ~ ^^'^ 

To be able to appreciate the relevance of our approximation, we define the distance dist 
between the predicted pole and the real pole as 

dist = ^ {MpA - Mp,,i,Y + {TpA - TpoieY • (4) 

This parameter dist helps us to see the rate of convergence of our sequence of 
approximants as is shown in Fig. 1 where N is the order of our approximation. In 
particular, the first prediction, using the P^, has an error below O.lMeV, the second 
one below 0.01 MeV, and so on. The reader should take into account that the plot 
is in logarithm scale. At that point, a word of caution is needed. To succeed on the 
prediction of the position of a resonant pole using the Montessus' theorem a crucial 
condition must be fulfilled. The resonant pole we are looking for must lie within the 
disk of applicability of the theorem, in our case, in the disk limited by both production 
thresholds. In that sense, the prediction of the position of the pole in a p-like model 
converge very fast, for a a-like model (which has a fatter resonance than a p-like model), 
the convergence is slower and for a ultra-fat-like model there are no convergence at all. 
This pattern of convergence is schematize in Fig. 2 where red dots represent a p-like 
model with a resonance at q^p = (0.77 — zO. 15/2)^GeV^, the blue squares a C7-Iike model 

with a resonance at q^p = (0.48 — i0.53/2)^GeV^ and the empty-green squares an ultra- 
fat model with a resonance at q^ = (2. 10 — il.05/2)^GeV^. All in all, we can conclude 
that the difference between the prediction of the q^ pole using the Pf and using the P^^^ 
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FIGURE 2. Rate of convergence of the Pade approximants sequence corresponding to the p-Hke 
model (red-circles), the a-like model (blue-filled squares) and the ultra-fat-like model (green-empty 
squares). 



log 



10 

-0.5 - 

-1.0 - 

-1.5 - 

-2.0 - 



dist 
IMeV 



2 3 4 5 



iV 



FIGURE 3. Rate of convergence of the fitted sequence Pade Approximants corresponding to a p-like 
model considered in the second case of application of our mehtod (see text for details). 



gives US an estimation of the systematic error of our method. 

As a second case of application of the Montessus' theorem, we generate a series of 
zero error data points for both modulus and phase-shift of our p-like model, Eq. (3). We 
then construct a generic Pade Approximant sequence which have several unknown 
parameters a„. To predict the position of the resonant pole we need to know the values 
of these a„ parameters. We fit then the modulus of our to the modulus of our p-like 
model and the phase-shift of our P^ to the phase-shift of our p-like model. Again, when 
the a„ parameters are known, we extract the position of the PA pole and we compare it 
with the real position through the distance dist. The rate of convergence of our new fitted 
sequence is shown in Fig. 3 and gives us again an estimation of the systematic error of 
our method in this second case. As expected, this convergence is slower compare with 
the previous genuine case. 



APPLICATION OF THE METHOD IN A PHYSICAL CASE 



We would like to use now our method to analyze the final compilation of ALEPH nn 
vector form factor data for the squared modulus |F;r;r(^^) [3], and the 7 = 7=1 nn 



TABLE 1. Comparison of different results for 
the determination of the mass and width of the 
p -meson. 



Ref. 


Mp(MeV) 


rp(MeV) 


[5] 
[6] 

[7] 
[8] 


762.5 ±2 
754 ±18 

763.0±0.2 
764.1 ±2.7+^;^' 


142 ±7 
148 ±20 

139.0±0.5 
148.2±1.9;^:^ 


this work 


763.7 ±1.2 


144 ±3 



scattering phase-shift 5^,:^:, identical to the TtTt vector form factor phase-shift in the elastic 
region 4m| <q^ < (if multipion channels are neglected), which conforms the range 
of applicability of our Pade Approximant sequence (more details can be found in Ref. 
[4]). For A'^ > 3 the fit already lies within the 68% confidence level (CL) and becomes 
statistically acceptable. At this point one needs to reach the typical fitting compromise. 
On one hand the experimental errors have an statistical origin and the contribution of 
this error increase as one considers higher order Pade Approximants , with larger 
number of parameters. On the other hand, the systematic theoretical error decreases as 
increases and the PA converges to the actual VFF. In the present work we have taken 
N = 6 ?& our best estimate as the new parameters of Pade Approximants with N >1 
turn out to be all compatible with zero, introducing no more information with respect 
to the previous Pf . Furthermore, the model studied before shows that the theoretical 
errors for mass and width results are smaller than 10~^ — 10~^MeV for A/^ > 6 (see 
Fig. 3), being negligible compared to the ^(1 MeV) experimental errors. This yields to 
the determination: 

Mp = 763.7 ± 1 .2MeV, Fp = 144 ± 3MeV , (5) 

which is found to be in good agreement with previous determinations shown in Table 
1 using more elaborated and complex procedures and with similar size of uncertainties . 

OUTLOOK 

A next step of this line of analysis would consist on incorporating another resonance 
in the problem where one should make use of the Montessus' theorem with a Pade 
Approximant sequence instead of the one (or, more general, P^ when a fixed number 
L of resonances are include in the problem). However, one should take into account that 
we are not demanding the function to be meromorphic. In that case, it is better to use 
the Pomerenke's theorem as was shown in [12], instead of the method developed here. 
On the other hand, since we are dealing with functions that present a branch cut and 
pole(s), one can make use of the so-called Multivalued Approximants, or Non-Rational 
Approximants, which generalize the usual Pade Approximants. Firstly considered by C. 
Pade himself, did not receive much attention until the works of Shafer [10] and Short 
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FIGURE 4. Rate of convergence of the QA Quadratic Approximant sequence corresponding to the 
p-hke model Eq. 3. 



[11]. The simplest of these approximants are the so-called Quadratic Approximants 
(QA), defined by 



(^) t + ^ W E + t = ^(^'■+'+'+') • 

n=0 n=0 n=0 

However, that can be extended to a given k order just following the definition: 



(6) 



A' 

n=0 n=0 



(7) 



where k = 2 would correspond to the QA, ^ = 3 to a Cubic Approximant and so on 
[ 1 1] . In our case, the order appearing in Fig. 4 is defined as N = l + k + j with N >2 
and / > 1 in order the QA to has a brunch cut. 

Such approximants possesses a build-in cut structure, they seem then ideally suited 
for problems related to multivalued functions such as our case. A preliminary study 
of the p-like model of Eq. (3) with this kind of Quadratic Approximants centered 
again at turns out to be very promising considering the their rate of convergence 
shown in Fig. 4 where we have selected the best approximant for each N among the 
N+{N-l) + {N-2)-{ hi possibilities. 

Finally, in several studies of S -matrix theory has been common to use, instead of 
the variable q^, a variable k related with the previous one through the relation k — 

(f- — Am^ (where m is the produced particle mass). In the fc-plane, the resonances 
are closer to the real axis than in the ^^-plane, so the extension of our method in this 
new variable suggests better pole position determination for those resonances with large 
width (such as the sigma meson). 



CONCLUSIONS 

We have develop a model-independent method for extracting resonance poles from 

physical amplitudes. The method is based on the well defined mathematical theory 
of Fade Approximants and makes use of the Montessus' theorem to systematize the 



algorithm of extracting the desired resonance pole. However, our method has a larger 
application since does not rely on a particular lagrangian or in a modelization on how to 
extrapolate from the data on the real energy axis into the complex plane (such as other 
methods discussed in Ref. [9]). In such a way, the analysis of other available form factors 
and phase-shifts will be presented elsewhere. 

In the particular case presented here we have analyzed the experimental Minkowskian 
;r;r-VFF and ;r;r-scattering data by a Pade Approximant sequence centered between 
the first and the second production thresholds (in such a way that we can unfold the 
second Riemann sheet of our amplitude). We have obtained a prediction of the p -meson 
pole position with compatible accuracy compared to other determinations of the same 
quantity (see Tab. 1) but with a simpler and systematic method. 
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